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Nomenclature 


AVT 

BART 

BL 

B.L. 

CAWAP 

CAWAPI 

CFD 

CFL3D 

Cf 

C P 

EFD 

ESP 

ET 

FC 

FS 

h 

HSCT 

HSR 

HUD 

ij,k 

ITAR 

LE 

Moo 

PSF 

PSP 

R n 


Applied Vehicle Technology 

Basic Aerodynamic Research Tunnel 

butt line on airplane, in., positive on right wing (See Fig. 1) 

Boundary layer 

Cranked Arrow Wing Aerodynamics Project 
CAWAP International 
Computational Fluid Dynamics 

name of structured grid flow solver developed at NASA Langley 

local skin friction coefficient 

Pressure coefficient 

Experimental Fluid Dynamics 

Electronically Scanning Pressure 

Exploratory Team 

Flight Condition 

fuselage station on airplane, in., positive aft (See Fig. 1) 

airplane altitude, ft 

High Speed Civil Transport 

High Speed Research 

Head-Up Display 

grid indices 

International Traffic in Arms Regulations 
leading edge 

free-stream Mach number 

Performance, Stability & Control and Fluid Physics 
Pressure Sensitive Paint 

Reynolds number, based on aircraft reference chord of 24.7 ft 
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RTB 

RTO 

V/Vre 

VFE-2 

VL 

VUB 

T 

WL 

x/c 

y 

+ 

y 

a 

P 

2y/bi 

2- D 

3- D 

Subscripts 

avg 

nom 


Research Technology Board 
Research and Technology Organization 

ratio of velocity magnitude in boundary layer to that at the Rake Extreme total-pressure tube 

Vortex Flow Experiment Number 2 

Virtual Laboratory 

Vrije Universiteit Brussel 

absolute temperature, °R 

waterline on airplane, in., positive up (See Fig. 1) 
fractional distance along the local chord, positive aft 
normal distance above the surface at a rake location, in. 

Reynolds number like term for flat-plate turbulent boundary layer 1 
angle of attack, deg 
angle of sideslip, deg 

fractional distance along the wing local semispan, positive toward the right wing tip 

Two-dimensional 

Three-dimensional 


average value 
nominal value 
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Introduction 


This paper provides a brief history of the F-16XL-1 aircraft, its role in the High Speed Research (HSR) 
program and how it was morphed into the Cranked Arrow Wing Aerodynamics Project (CAWAP). Various flight, 
wind-tunnel and Computational Fluid Dynamics (CFD) data sets were generated during the CAWAP 2 . These are 
described and sample data comparisons given. This is followed by a description of how the project became 
internationalized and is concluded by an introduction to the results of a four year CFD predictive study of data 
collected at flight conditions. 

CAWAP History 
F-16XL airplane 

The F-16XL-1 airplane is a single-place, fighter-type prototype aircraft developed by the General 
Dynamics Corporation-Ft.Worth Division (now Lockheed Martin Aeronautics Company - Fort Worth) by stretching 
the fuselage of a Full-Scale Development F-16A and adding a cranked-arrow wing, a modified fuel system, and a 
modified flight control system. There were two F-16XL aircraft built, the -1, which was used in the CAWAP, and a 
two-place version, the -2. Consequently, the terms F-16XL and F-16XL-1 are used interchangeably in this paper. 
Both aircraft had scheduled leading-edge flaps, elevons, and ailerons on the wing for control. The technical 
specifications for the airplane are given in Table 1, which is reconstructed from reference 2. Details on the 
construction of the aircraft and its intended missions are given in references 3-5. 
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Table 1. F-16XL-1 Airplane Specifications 


Feature 

Value 

Wing Span 

32.4 ft 

Height 

17.606 ft 

Length 

54.155 ft 

Reference Chord 

24.7 ft 

Theoretical Root Chord 

41.75 ft 

Wing Area 

646.37 ft 2 

Reference Wing Area 

0 \ 

O 

O 

Reference Aspect Ratio 

1.75 

Typical Takeoff Weight 

35,000 lb 

Engine; Max Thrust 

Pratt & Whitney F100-PW-200; 23,830 lb 


The design of the cranked-arrow wing was a cooperative effort of the General Dynamics Corporation and 
the NASA Langley Research Center. The new wing, common for both the -1 and -2 versions, was designed to 
provide the F-16 aircraft with improved supersonic performance while maintaining comparable transonic 
performance to that provided by the current F-16 design. As shown in Figs. 1 and 2, the resultant design had a 
leading-edge (LE) sweep angle of 70° inboard and 50° outboard of the crank. At the juncture of the wing leading 
edge with the fuselage, an “S-blend curve” was placed in the leading edge to alleviate a pitch instability that was 
found to occur at high angles of attack in wind-tunnel tests. All flight test data reported in reference 2 were 
collected with the air dams - upper-surface fences mounted near the wing leading-edge crank - and wing-tip 
missiles installed, as shown in the figures. 
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Figure 1. Three-view drawing of F-16XL-1 airplane. 



Figure 2. F-16XL-1 aircraft with missiles, tufts, modified flow-visualization paint scheme, and video targets 

(NASA Photo EC96-43508-2). 
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HSR program overview 

This program was started in 1 990 to develop the technologies that would result in a supersonic passenger 
jet that would fly 300 passengers at more than twice the speed of sound while maintaining comparable ticket prices 
to current subsonic transports. As envisioned by the government and industrial partners, the high speed civil 
transport (HSCT) would cross the Atlantic or Pacific oceans in half the time of modem subsonic jets using new 
technologies for airframe manufacturing, propulsion systems, aerodynamics, and reduced environmental impacts. 
By 1995, based on several industry design concepts, computer modeling, and wind tunnel tests, a Technology 
Concept Airplane was selected as a common reference point in furthering the technology development process. This 
single concept was to have improved aerodynamic performance and operational characteristics while also meeting 
environmental goals for emissions and noise pollution. The HSCT concept is depicted in Fig. 3. 



Figure 1. High-Speed Civil Transport concept. 


Due to economic constraints, the program was phased out in fiscal year 1999. Despite the cancellation, the 
program still managed to meet or exceed many of the original HSR program goals. New high-temperature 
composite materials and structural concepts were developed to keep weight and manufacturing costs down while 
maintaining the strength and durability that would be required for the aircraft. NASA engineers developed new 
vision systems for the pilots that maintained safety and performance capabilities while eliminating the need for a 
drooped nose similar to the Concorde. Another critical element to the program was the development of a propulsion 
system that would not harm the environment both in atmospheric effects as well as mitigating the noise, while 

providing the performance and durability required to keep the aircraft economically viable. 
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F-16XL flight overview and planning for HSR 


Early in the development of the HSR program, it was determined that improved aerodynamic performance, 
while reducing the noise at high-lift conditions during take-off and landing, was a critical element to the program. 
The availability of the F-16XL aircraft with its cranked wing configuration, which was similar to the HSCT concept, 
provided a unique opportunity for CFD correlation and code validation with flight and wind-tunnel data. The 
project was divided into three phases based on the required geometry changes to the basic aircraft. The phases are 
depicted schematically in Fig. 4 relative to the aircraft modification required. Phase 1 was the basic airplane with no 
modifications. This phase would serve as a baseline for the future modifications. Phase 2 required the removal of 
the “S” curve in the apex region of the wing. This second configuration was more representative of the HSCT 
cranked wing. The third and final phase would incorporate a high-lift device along the entire leading edge of the 
aircraft. The exact configuration of the high-lift device would be determined from wind-tunnel experiments and 
CFD predictions. The high-lift configuration would be representative of the HSCT in a take-off or landing 
configuration. 



Figure 2. F-16XL-1 research aircraft modifications. 

The primary objective of the flight test project was to verify the performance of high-lift concepts while 
ensuring compliance with community noise standards. The first step was to establish a ground to flight correlation 
for the cranked-arrow wing planform. All three configurations as described would be flight-tested and the data 
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would be used as a calibration of the design analysis tools as well as the noise prediction codes. In addition, 
advanced operating procedures for take-off and landing would be evaluated. The final objective was to assess 
integration and “real-world” operation of high-lift devices. A schematic of the ground to flight correlation is shown 
in Fig. 5. 
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Figure 3. Schematic of ground to flight correlation for assessing the impact of aerodynamic and operational 

concepts on noise profiles. 


CAWAP 

During the Phase 1 experiments on the baseline F-16XL aircraft in the spring of 1994, the remainder of the 
project as planned was cancelled due to funding limitations. Sufficient funding was provided to complete the first 
phase of the project, with a slight change to the objectives. This was the start of the CAWAP as it is known today. 
The revised objectives were to document upper-surface flow physics at high-lift and transonic test conditions and to 
characterize the stability and control of the aircraft. The original intent of a flight, wind-tunnel, and CFD correlation 
experiment would be maintained, albeit for the baseline F-16XL configuration only. Table 2 illustrates the 
extensive set of planned comparisons between flight, wind-tunnel and CFD, and Table 3 provides the actual set of 
comparisons made. (These two tables are reconstructed from Ref. 2 with modifications shown in italics from Ref. 2. 
The modifications are for additional comparisons that were overlooked at the time of the publication.) 
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Table 2. Planned F-16XL CAWAP Data Comparisons 


Item 

Data comparison 

Data type/instrument 

Data source 

1 

On and off surface flow 

Tufts 

Static pressures in a row 
Vapor screen 
Particle traces 

Flight 

Flight 

Flight 

CFD 

2 

Surface flow and pressure 

Oil flow 

Pressure sensitive paint (PSP) 

Wind tunnel 
Wind tunnel 

3 

Surface pressure 

Pressure sensitive paint 
Calculations 

Wind tunnel 
CFD 

4 

Vortex core location 

Vapor screen 
Particle traces 

Wind tunnel 
CFD 

5 

Vortex core location 

Vapor screen 
Particle traces 

Flight, wind tunnel 
CFD 

6 

Surface flow 

Oil flow 
Tufts 

Liquid Crystals 

Flight 

Flight 

Flight 

7 

Surface flow 

Oil flow 
Tufts 

Particle traces 

Flight 

Flight 

CFD 

8 

Surface flow 

Oil flow 
Particle traces 

Flight, wind tunnel 
CFD 

9 

Surface pressure 

Pressure sensitive paint 
Particle traces 

Wind tunnel 
Wind tunnel 

10 

Off surface flow 

Vapor screen 
Five-hole probe 

Flight, wind tunnel 
Wind tunnel 

11 

On and off surface flow 

Oil flow 
Vapor screen 
Particle traces 

Flight 

Flight 

CFD 

12 

On and off surface flow 

Oil flow 
Vapor screen 
Five-hole probe 

Flight, wind tunnel 
Flight, wind tunnel 
Wind tunnel 

13 

Surface flow 

Oil flow 
Tufts 

Liquid Crystals 

Propylene-glycol-methyl-ether traces 

Flight 

Flight 

Flight 

Flight 

14 

Surface flow 

Oil flow 
Tufts 

Particle traces 

Propylene-glycol-methyl-ether traces 

Flight 

CFD 

Flight 

Flight 

15 

Vortex core location 

Vapor screen 
Dye traces 
Particle traces 

Wind tunnel 
Water tunnel 
CFD 

16 

Surface pressure contours 

Static pressures 
Calculations 

Flight, wind tunnel 
CFD 

17 

Surface pressure 

Pressure sensitive paint 
Static pressure in a row 

Flight 

Flight, wind tunnel 

18 

Surface pressure 

Pressure sensitive paint 
Static pressures in a row 
Static-pressure surface 

Flight 

Flight, CFD 
CFD 

19 

On and off surface flow 
contours 

Tufts 

Particle traces 
Stagnation pressure 

Flight 

CFD 

CFD 
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Table 2. Concluded. 


Item 

Data comparison 

Data type/instrument 

Data source 

20 

Boundary-layer profile 

Rake 

Velocities 

Flight 

CFD 

21 

Skin-friction distribution 

Modified Preston tube 
Equation 

Flight 

CFD 

22 

Leading-edge boundary layer 

Hot-fdm gauges 

Flight 


Table 3. Actual F-16XL CAWAP Data Comparisons 


Item* 

Data comparison 

Data type/instrument 

Data source 

3 

Surface pressure contours 

Pressure sensitive paint 
Calculations 

Wind tunnel 
CFD 

6 

Surface flow 

Oil flow 
Tufts 

Flight 

Flight 



Liquid Crystals 

Flight 

12a 

On and off surface flow 

Oil flow 
Vapor screen 

Wind tunnel 
Wind tunnel 

16a 

Surface pressure contours 

Static pressures 
Calculations 

Flight 

CFD 

17a 

Surface pressure 

Static pressure in a row 

Flight, wind tunnel 

18a 

Surface pressure 

Static pressures in a row 
Static-pressure surface 

Flight, CFD 
CFD 

19 

On and off surface flow 

Tufts 

Particle traces 

Flight 

CFD 



Stagnation pressure contours 

CFD 

20 

Boundary-layer profile 

Rake 

Velocities 

Flight 

CFD 

21 

Skin-friction distribution 

Modified Preston tube 
Equation 

Flight 

CFD 

22 

Leading-edge boundary layer 

Hot-film gauges 

Flight 


*The “a” in the item number signifies a reduction in the actual versus planned number of items being compared. 


Note that in Table 2 an attempt was to be made to perform off-surface laser- vapor-screen data using seeded 
material, and to develop surface streaklines using propylene-glycol-methyl-ether. However, these two types of data 
were not obtained due to the lack of funding and higher priority datasets that needed to be collected. Moreover, 
consideration was given to perform PSP in flight; however, it quickly became apparent that the timing and cost of 
the technological development was beyond the scope of this project. The concept was to use the laser light-sheet in 
combination with a surface coating to obtain the pressure data, but it only reached the idea stage in a discussion 
between NASA Langley researchers and key personnel from the airframe industry. In addition, since the laser 
system was not funded to completion, this idea had to be abandoned. 
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Data available from CAWAP 


Seven different types of flight data were collected, as per Table 3, and four are shown schematically in 
Fig. 6. Three were pressure based - surface static pressures, boundary-layer rakes, and modified Preston tubes 6 ; 
three were video-recording based - surface tufts, surface oil and surface liquid-crystals; and one was hot-film data. 
The pressure- and surface-flow-data are used for the purpose of establishing the effects of variation in Mach number 
on the local flow. These data serve as the basis for comparison with other data sets. The hot-film data are used to 
establish whether boundary-layer transition occurs and under what test conditions. Geometry data of the airplane 
upper surface was also obtained, using photogrammetry, and compared with the numerical surface description (See 
Ref. 2). 



Figure 4. Layout of pressure & hot-film instrumentation suite on the F-16XL-1 aircraft. 


Pressure suite 

Figures 6 and 7 detail the complete pressure instrumentation system layout on the aircraft, including the 
distribution of the static ports by type, belt or flush, and boundary-layer rakes or modified Preston tubes. These 
static ports are connected to internally mounted Electronically Scanning Pressure (ESP) modules through tubes. 

12 

American Institute of Aeronautics and Astronautics 



Each pressure-tube in the belt was used to measure two separate values of pressure. This was accomplished by 
sealing each tube about halfway along its length, thereby making provision for one forward and one aft port. The 
actual number of static ports associated with each belt are displayed in Fig. 7 imbedded along the belt nominal Butt 
Line (BL). 



Figure 5. Details of complete-pressure-instrumentation suite and layout on the airplane. 

The right-wing surface pressures - mostly upper surface - were measured using 337 static ports, both flush 
(LE region) and in streamwise belts, through eleven 32-port ESP transducers - also called modules. Of these 337 
ports, only 326 proved to be reliable and the distribution was 280 on the upper surface and 46 on the lower. The 
ports were arranged so that there would be a sufficient number at a given BL or Fuselage Station (FS) for cross 
plotting, as well as for covering other regions of special interest; i.e., the apex and ahead of/and behind the hinge- 
lines of the trailing-edge control-surfaces. 

Boundary-layer measurements were made using two, two-inch high, rakes at a time at four different 
positions on the left wing and the most inboard one was always used as a control. Each rake used 16 active tubes - 
15 total pressures and one static pressure - of the 23 available. These two rakes were connected to one 32-port ESP 
module located inside that wing. When mounted on the aircraft, each rake was oriented into the local flow at an 
average angle over its height based on initial CFD predictions from the CFL3D code. 7 ' 8 The flow conditions were 
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for the complete aircraft (half-aircraft modeled with symmetry assumed) at a = 13°, M 00 = 0.29 and R n = 46.1xl0 6 ; 
i.e. Flight Condition (FC) 7. (See Tables 4 and 5 for a listing of the FCs of interest in the current study and the 
associated flight and engine parameters, including those for FC7.) 

Figure 8 shows the four boundary-layer locations chosen - one well inboard of the shed vortex systems 
(#3), one underneath the primary vortex (#4), and two associated with the secondary vortex. One of the latter is 
located underneath that vortex (#5) and the other (#7) at its separation point; all are at a nominal position of FS 295. 
The average of the local flow at- and slightly off-the- surface were used to establish the rake orientation angles for 
FC7. This figure also shows the relative locations of the modified Preston tubes. They were to be located at the 
same fuselage station as the boundary-layer rakes, but a more aft position for the tubes was necessitated due to 
easier aircraft installation and to avoid the flow off a step in the leading-edge region. 

The 16 modified Preston tubes (See Figs. 6 and 8) - the modification to each Preston tube is the integration 
of a static pressure port with the total pressure tube - are used to determine local skin friction across the left wing 
near FS 330. These 32 pressures use the same ESP as the rakes but not on the same flight. The tubes were aligned 
with the local flow using the same initial CFL3D solution at FC7. The equation used to generate the Experimental 
Fluid Dynamics (EFD) c f values comes from reference 5 and relates, through a process of calibration, the pressure 
change between the total and static tubes to the local skin friction. 


Figure 6. General arrangement of rake and modified Preston tube relative locations on F-16XL-1 left wing; 
Pressure instruments oriented for Moo= 0.291, a= 13°, and R n = 46.1xl0 6 


Boundary layer rake 



#3 i— Extension of hot -fi lin gauge 
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Video suite 


Video data was recorded with up to six external cameras; two mounted atop the vertical tail, one on either 
side of the fuselage behind the canopy, and one in the nose of each dummy missile. An internally mounted HUD 
camera was also used on occasion. Figure 9 shows the camera locations on the aircraft. The time was added to each 
image by a time-code inserter so that the images could be compared to form a composite and so that the flight test 
conditions could be established. Images of interest were digitized in a 512 by 480 pixel format for further 
processing in order that quantifiable video data be developed. In addition to the images, the other input quantities 
needed for the processes are the video targets and the position and calibration characteristics of each camera/lens 
combination. 



Figure 7. Locations of the video suite on the F-16XL-1 airplane. 


Wind Tunnel 

The second major source of data was from wind-tunnel tests with a variety of model scales. Some of these 
wind tunnel tests were conducted on the F-16XL configuration prior to the HSR program and others were done in 
direct support. These tests encompassed a variety of configurations at various Mach and Reynolds numbers. Since 
some of the tests were conducted with the air dams not installed, they were not applicable for comparison with flight 
data. The following three wind tunnel tests were primarily used for the flight and CFD comparisons. 


15 

American Institute of Aeronautics and Astronautics 


1. A test on a 0.1 1 -scaled model was conducted in the NASA Ames 1 1-Foot Tu nn el prior to the start of the 
flight project. The test is described in reference 9, while the actual data is tabulated in reference 10. This model was 
specifically built to estimate the airloads for the airplane from M=0.60 to 2.0 using 190 pressure ports distributed on 
the left-wing upper surface and the right-wing lower surface. The Reynolds numbers tested were 2.3x1 0 6 and 
2.75xl0 6 , while the angle of attack varied from -1.94° to 28.75°. The ports were located in streamwise rows that did 
not nominally match the airplane. Because of data release restrictions, no direct comparisons were made with flight 
or CFD results. However, un-scaled transonic pressure coefficient contours were used for comparison purposes in 
reference 2. 

2. The test for the 0.18-scaled model was conducted in the NASA Langley 30-by 60-Foot Tunnel at a Mach 
number of less than 0.08 and a Reynolds number of 2.1xl0 6 . The angle of attack varied from -5° to 30° and the 
angle of sideslip varied from -20° to +20°. In addition to force/moment and pressure data, some vapor screen 
images were captured. A set of pressure results was published 2 and the force/moment results were published in an 
earlier paper. 11 The basic model has 30 right-wing, flush upper surface ports located to yield pressures for both 
streamwise and span wise rows. The ports are nominally duplicated on the airplane for comparison. Only two runs 
were utilized from this test for comparison. One was a force/moment and the other a pressure, both at 0° sideslip. 

3. The third wind-tunnel test used for comparison was a 0.04-scaled model in the NASA Langley Basic 
Aerodynamic Research Tunnel (BART). The test conditions included Mach numbers less than or equal to 0.165, 
Reynolds numbers less than 1.12xl0 6 , and an angle of attack variation from 5° to 20°. This model had 82 pressure 
ports divided between the right upper wing surface and the left lower wing surface. The pressure ports correlated 
with the 0.18-scaled model as well as being nominally duplicated on the airplane. Selected results from the BART 
test appear in reference 12, where they have been compared with Euler code predictions. 

CFD modeling 

The last major source of data was that generated by CFD 2 . The flow solver CFL3D was run in the Navier- 

Stokes mode with a turbulent boundary-layer employing the Baldwin-Lomax with the Degani-Schiff turbulence 

model (in the j-k directions) on a multiblock, patched grid over a variety of wind-tunnel and flight test conditions 

(See Refs. 7, 8). Two separate grids were used to model half the aircraft configuration (with undeflected control 

surfaces) and external flow field. The initial grid had 36 blocks and was used with version 3 of the flow solver to 
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produce the initial results upon which the locations and orientations of the surface instrumentation were set. The 
current grid had only 30 blocks and was used with version 5 of the flow solver to obtain the comparative solutions 
reported herein. The current grid was needed for two reasons; namely, to have the grid more closely conform to the 
actual fuselage and wing geometries, and to improve the grid layout on the wing and fuselage surfaces. The missile- 
and missile-rail-grids were effectively the same with either grid. For the current grid, the inner region of the aircraft 
was modeled by 16 blocks, the outer region by 14 blocks, and all 30 blocks are shown schematically in Fig. 10. The 
boundary conditions were symmetry, solid wall for the outer mold lines, flow into the duct-inlet but the exhaust-face 
was faired over, and Riemann-type conditions at the far-field boundaries. A total of 1,372,096 cells (1,707,1 17 node 
points) were used to obtain solutions at specified test conditions (M*,, a, grid R n , T, etc.). To maximize computer 
resource allotments, the minimum number of cells was used. The resulting grid spacing normal to the numerical 
surface led to a value of y + of 2 at wind-tunnel R n , whereas at flight R n the average value was y + of 82. In an effort 
to compensate for the insufficient grid spacing at flight conditions, the “wall function” option was used to augment 
the turbulence model in CFL3D. The “wall function” is defined as that boundary-layer growth rate expected from a 
turbulent mean flow near the wall 13 . 



Fuselage 
Inboard upper wing 



(a) Inner. 
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(b) Outer 

Figure 8. CFD block structure layout for F-16XL-1. 


Database of results 

A database was set up, as described in Appendix B of reference 2, to facilitate data comparison between the 
various CAWAP sources; i.e. flight, wind-tunnel and CFD. The data types used for comparison were local surface 
pressures between the three sources, and boundary-layer and skin friction between the flight data and CFD 
predictions. Moreover, the database stores administrative information about the tests (meta-data) and the full, mass 
storage, path name of the resulting data files. For the EFD measurements, the data stored are files of pressure, 
force/moment, still photographs, and two-dimensional images digitized from videotape (both instantaneous and time 
averaged); for CFD predictions, the grid and solution files are stored as well as C p data at selected FSs and BLs. 
The F-16XL-1 database also supports the viewing of 3-D renderings of the 2-D flight images through computer 
software tools. Queries for selected pressure data and 2-D flight image data were available to internet users by 
completing an appropriate on-line form on this unclassified/unlimited server. The web service was functional at the 
publication time of reference 2 and shortly thereafter, but has since been discontinued due to changes in policies and 
available resources. 
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Sample CAWAP comparisons 


Selected figure samples from the CAWAP report 2 are repeated as Figs. 11 to 14. Following each set of 
figures is an associated commentary provided for completeness. Note that these figures may contain a legend 
signifying which flight and test point (run number) served as a source for the flight data; for example, 
Fltl44.Run3b, up. These dataset associations are not the same as the FC but each one shown has a FC equivalent, 
as per the figure caption. 



(c) BL 80. 


(d) BL 95. 
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(g) BL 184.5. 


(h) FS 185. 
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(k) FS 375. 


(1) FS 407.5. 



(m) FS 450. 


(n) FS 492.5. 



(o) Upper surface Cp distribution. 

Figure 9. CFD (CFL3D) and measured Cp at FC46 (M. = 0.527, a = 10.4°, R n = 46.9 X 10 6 ) from Ref. 2. 
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Figures 11(a) to ll(n) show good overall agreement between the measured and predicted results at this 
flight condition. There are, however, three notable exceptions. The first is that the primary vortex effect on the 
suction peak at BL 55, the most inboard location, is under-predicted. The second is at BL 153.5 where the measured 
data has its primary- vortex suction-peak more forward than predicted and is followed by a more rapid compression 
downstream. Lastly, is the unusual and uncharacteristic variation near 2y/bi=0.6 for FS 185. The associated ports 
for the latter are located in and towards the aft-end portion of the S-blend curve part of the airplane, a region where 
the flow is very sensitive. Even with this unusual behavior of C p , these results are retained because they do bracket 
the predicted data at this FC. 

Figure ll(o) provides an overview of the upper surface C p comparison where the black dots represent the 
location of the flight ports and the associated color of the surrounding bubble indicate the measured value. Where 
the color of the bubble is the same as the background CFD solution, only the black dot is noted. Generally good 
overall agreement is seen except in the S-blend region where the measured pressures are more negative, but 
elsewhere the differences are such that the measured values are more positive. 

Figure 12 shows the fusion (overlaying) of surface tuft images from the three left-side cameras (tail, 
missile-nose and fuselage-apex) with CFD surface streamlines and vortex-core representation for FC46. [The white 
circles are video targets used for camera registration.] In particular, Fig. 12(a) presents the camera combination 
obtained by using the video targets with camera location and orientation angles. Figure 12 (b) shows that the CFD 
surface streamlines compare well with the flow depicted by the surface tufts. Figure 12 (c) shows the addition of the 
stagnation pressure iso-surfaces - at a value of 0.78 representing the locations and extent of the various airplane 
vortex systems - with transparency. As expected from the results of the surface comparison, the vortex system is 
well located with respect to the flight tufts 2 . 


22 

American Institute of Aeronautics and Astronautics 





(a) Tuft images projected from three cameras 
onto aircraft grid. 


(b) Combination of tuft images 
and CFD surface streamlines. 



(c) Combination of tuft images, streamlines, and vortex systems. 


Figure 10. Flight tuft data from three cameras on F-16XL-1 airplane and CFD solution at FC46 (M* = 0.527, 

a = 10.4°, Rn= 46.9xl0 6 ) from Ref. 2. 


Figure 13 shows the boundary-layer profiles of measured and predicted data. The profile for B.L. rake #3 
is seen in Fig. 13(a) to be that of a classical, streamwise flow and this location was measured on all such flights to be 
used as a control. This figure shows good data repeatability and that the predictions agree well with them away 
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from the wall. The data presented in Fig. 13(b) for B.L. rake #4, underneath the primary vortex, shows that the 
predicted trends are seen in the measured data, including the y location for the onset of ‘jet-like’ flow followed by a 
reduction from there to the rake extreme. Figures 13(c) and 13(d) are associated with the secondary vortex, one 
underneath the core (#5) and one along the secondary separation line (#7). Two points are highlighted with respect 
to these two rakes data sets: the predicted values do not match the measured ones but do show how the anticipated 
flows would impact the boundary layer profiles; and for all practical purposes the two measured data sets are the 
same with a profile that does not asymptote to the boundary-layer rake extreme value. 
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(d) Rake 5: FS 294.59, BL -96.06. 


Figure 11. Predicted and measured velocity profiles for boundary-layer rakes on F-16XL-1 airplane for FC7 
(M.,, = 0.29, a = 11.89°, p nom = 0°, h = 5000 ft, R„ = 44.4 x 10 6 ) from Ref. 2. 


24 

American Institute of Aeronautics and Astronautics 






Figure 14 shows the predicted and measured local skin-friction coefficient at FS 330 for FC7 and FC19, 
respectively, as these FCs were quite close. High skin-friction peaks are an indication that a vortex system is present 
above the surface. Both sets of data show at least two such peaks. The predicted values for the primary (inboard) 
vortex have a different trend and reach a lower value than the measured data, and the vortex location is predicted 
more inboard than that measured - ~BL -90; however, the solution does a somewhat better job in estimating the 
peak value and location of the secondary vortex, albeit with two unexplained oscillations which trail toward the LE 
(BL~ -116). 



Figure 12. Predicted and measured c f on F-16XL-1 airplane at FS 330 for FC19/FC7 (M avg = 0.33, 
a avg = 11.9°, p nom = 0°, R„, avg = 45.6 X 10 6 ) from Ref.2. 

CAWAPI 

Organization 

At the Spring 2000 Research and Technology Organization/Applied Vehicle Technology (RTO/AVT) 
Symposium week held in Braunschweig, Germany, this international group was asked whether there was interest in 
predicting the flight results around the F-16XL aircraft - as the report 2 documenting the work was in the process of 
being readied for publication. As a consequence, the Performance, Stability & Control and Fluid Physics (PSF) 


Technical Committee established an Exploratory Team (ET) to assess interest among the member nations on this 
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and related topics. During the ensuing year, a team was constituted and, at their first meeting in Loen, Norway, in 
the spring of 2001, they examined a variety of vortical flow topics that could lead to task groups or symposia. A 
total of five topics were offered for consideration. Some went forward individually, while others were combined in 
order to gain AVT Panel recommendation for acceptance and Research and Technology Board (RTB) approval. 

In particular, two topics presented at the Loen Symposium on “Vortex Flows and High Angle of Attack”, 
during the 2001 spring RTO meeting week, were put forward separately for further study. They were the CAWAP 14 
and the Vortex Flow Experiment - Number 2 (VFE-2) 15 - both of which dealt with vortical flows around slender 
wings. In the fall of 2002, these two topics were merged into a single proposal that the AVT Panel could 
recommend for approval by the RTB. Success was achieved in the spring of 2003 when the RTB granted approval 
for AVT-1 13 and authorized its first meeting to commence during the upcoming RTO/ AVT Symposium week. The 
title of the task group is “Understanding and Modeling Vortical Flows to Improve the Technology Readiness Level 
for Military Aircraft”. Even though integrated, each topic had its own facet lead, and these leaders became the co- 
chairmen of the task group. The title for the expansion of CAWAP activity to the International aerodynamics 
community within the RTO was denoted as CAW API. This name was crafted in advance in anticipation of approval 
and involved using a Virtual Laboratory (VL) at NASA Langley for facilitating secure data storage and transmission 
(See Refs. 16, 17) 

Features 

The features of the CAW API are contained in the foundational document for the task group known as the 
“Terms of Reference” - an internal RTO/AVT panel document. In particular, there are basically three objectives to 
be performed under this task among the participants. They are detailed as follows: 

1. Assessing various CFD codes against F-16XL-1 flight, and perhaps wind-tunnel, data sets in order to 
increase the Technology Readiness Level of the respective codes to a value of 5 (“Component and/or 
breadboard verification in a relevant environment”). 

2. Developing best practices for each code based on the data sets. 

3. Incorporating appropriate or upgraded turbulence models into the respective codes to provide for improved 
agreement. 
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These were to be accomplished by having each of the participating groups be responsible for certain 


aspects of the work. For example, each participating group will use their best efforts, consistent with program 
priorities and funding, to perform the agreed upon detailed tasks, and to be responsible for providing its own 
resources for the completion of this task. In particular, NASA agreed to do the following: 

1. Supply export-controlled geometry in various formats (iges, structured grid and unstructured grid) of the 
F-16XL-1 aircraft to participating partners once formal Memorandum-of- Agreements are in place. 

2. Make available flight pressures, images, skin-friction and boundary layer measurements to the team. 

3. Supply data format to and coordinate database services needed by the group. 

4. Coordinate the various efforts and arrange for meetings each 6 months in conjunction with RTO 
Symposium. 

5. The completion date of this task should be a maximum of 3 years after all approvals are granted. The 
initial ending date was set for December 2005, but this was later extended to December 2007, due to a 
variety of factors. 

The fact that the aircraft geometry was restricted by ‘International Traffic in Arms Regulations’ (ITAR) 
complicated the process of NASA fulfilling its obligations until it was realized that a Virtual Laboratory, housed in 
an electronically secure DeMilitarized Zone, was a solution for geometry and grid transfers among participants. 
Details about this have been documented 16,17 . Accessibility to both old and new data was provided through the VL 
to the participants. After the conclusion of the task group, it was anticipated that portions of the database would be 
made widely available via the VL, but that is highly unlikely given the current environment. 


The CAW API facet work was facilitated by having well-known and highly respected 
organizations/researchers as members of this international effort working under the RTO “umbrella”. These 
organizations appear in the following list grouped according to their basic function. They include four airframe 
companies — Boeing-St.Louis Phantom Works (USA), EADS-Munich (Germany), Lockheed Martin Aeronautics 
Company - Fort Worth (USA), and Turkish Aircraft Industries (Turkey); two government research laboratories — 
NASA Langley (USA) and NLR (The Netherlands); and six university led efforts — KTH/FOI (Sweden), United 
States Air Force Academy (USA), University of Glasgow/Liverpool University (United Kingdom), University of 
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Tennessee-Chattanooga/SimCenter (USA), University of Wyoming (USA), and Vrije Universiteit Brussel 
(VUB)/NUMECA (Belgium). Unfortunately, not all of these were able to continue to the end of this facet, and 
some optioned to continue in the VFE-2 facet instead. 


Presentations 

The original CAWAP solutions were obtained using a structured grid solver, based on a documented iges 
file with refinements, compared with measured flight data, and reported in reference 2 for the F-16XL-1 aircraft. It 
was anticipated that the new solutions would employ both structured and unstructured grids. Rather than just use the 
same iges file as before, it was decided to reinvestigate the geometry, find the best iges file available, and make 
certain that it was equally suitable for both solver types. The new iges file is only slightly different from the 
previous one, but it satisfied both ‘grid’ communities. The process of obtaining grids for both structured, as well as 
unstructured, solvers from this iges file is detailed in reference 18 and Fig. 15 is a sample from that paper. 




Figure 13. Abstraction of the surface geometry (left) and projected abstraction (right) for the F-16XL half- 
span model from Boelens,et.al. 18 . [Structured (pictured) and unstructured grid development from the iges 

file is discussed in this paper.] 

When the CAWAPI computational effort commenced, the majority of those performing computations 

planned to do so using structured grid solvers; however, in the ensuing years that trend has reversed and now the 

majority are using unstructured grid solvers. A partial reason for this movement has been the improvements made in 

these solvers that include the potential and actualization of breaking a solution down into components for parallel 

processing, and the potential for automatic grid generation. Both of these can lead to significant time reductions 

from geometry specification to solver results. In the listing of papers that follow, authors have utilized their solvers 
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- structured 19 " 21 or unstructured 22 ' 28 - in order to determine the best practices for each code and type of flow 
combination and then made comparisons with each other’s predictions or flight data as they deemed appropriate. 
Each of these papers has its own focus, but also contains common comparisons. For example, reference 19 
emphasizes surface streamlines and understanding the relation with the off-surface flow, as shown by the sample 
Fig. 16. A variety of aircraft modeling studies were performed in reference 20 with the results for the complete 
aircraft shown in Fig. 17. References 21, 24 and 28 emphasize the effects of various turbulence models and are 
shown by sample Figs. 18, 19, 20, respectively. Figure 20 28 also shows, for comparison, the flight data represented 
by small circles surrounding the pressure port locations, shown as black dots, and colored by the CFD C p scale. 
Reference 22 emphasizes solution adaptive grids, as shown by the sample Fig. 21. References 23 and 25 highlight 
time-accurate and time-dependent solutions, as well as time averaged, and are illustrated in Figs. 22 and 23, 
respectively. Figure 24 26 shows Mach 1 iso-surface colored by total pressure at FC70. Figure 25 27 shows 
streamlines and the x-vorticity component for FC50 where p = +5.31°. Lastly, reference 29 provides a summary of 
the “Lessons Learned” as a result of this computational effort, including its impact on the Technology Readiness 
Level of the examined current solvers. 

Please note that all of the above referenced papers have been summarized into the following five articles in 
this “Prediction of F-16XL Flight Flow Physics” special section. These articles follow the same general 
arrangement as the conference papers in that the next paper in this sequence is the geometry/grid discussion (ref. 
30), followed by solutions for the structured grid (ref. 31), unstructured standard grid (ref. 32), unstructured tailored 
grid (ref. 33), and a concluding article on the lessons learned (ref. 34). Each of these articles has a lead author who 
planned the scope and coordinated the efforts of the co-authors. 
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Figure 14. Surface streamlines with interpretation for FC7 with ENSOLV solver from Boelens, et.al. 19 . 


An initial group of four flight conditions with either vortex-dominated or transonic flows were adopted by 
this facet for prediction and comparison. These were later expanded to seven and included two sideslip conditions 
(See Table 4). Associated with each was a set of pressure/temperature/Mach number values for a generic engine 
(See Table 5). 



Figure 15. Upper and lower surface C p results for FC7 with PMB solver from Badcock 20 . 
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Figure 16. Turbulence modeling effect on C p for FC46 with PAB3D solver from Elmiligue, et.al. 21 . 


T — r — i — r" 


0.008 


0.006 - 


Jf 0.004 


0.002 - 


m 

:'A 

t C ■ 


r ; 1 T f 1 r f 

□ Flt131.Run14B. FC19 

CFL3D@FC07 

S-A 

EARSM 
CC -EARSM 
• “ ■ EARSM Hasten 
DRSM 



,0 fe 


* 100,0 


- 75,0 
BL in 


- 50,0 


- 25,0 


Figure 17. Turbulence modeling effect on C p for FC19 with Edge solver from Goertz, et.al. 24 . 
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FC46 Non-linear turbulence 


Y 



Figure 18. Turbulence modeling effect on C p for FC46 with USM3D solver from Lamar, et.al. 28 . (Flight data 
represented by small circles surrounding the pressure port locations, shown as black dots, are colored by the 
CFD C p scale. Agreement between the two sets is good when no discernable color difference is noted.) 



Figure 19. Solution adaptive grid effect for FC46 with TAU solver from Fritz 22 . 
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Figure 20. Time accurate solution for FC19 with Cobalt solver from Morton, et.al. 23 . 
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Figure 21. Time dependent solution for FC7 with BCFD solver from Michal, et.al. 25 . 
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Figure 22. Mach 1 iso-surface colored by total pressure at FC70 with Falcon v4 solver from Davis, et.al. 26 . 



Figure 23. Pressure-colored streamlines and x-vorticity component for FC50 with TENASI solver from 

Karman, et.al. 27 . 
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Table 4. Seven Flight Conditions to be Examined. 


Flight Condition 

Actual Mach 
Number 

Actual a, 
degs 

Actual p, 
degs 

Actual Reynolds 
Number 

FC7 

0.304 

11.89 

-0.133 

44.4x1 0 6 

FC19 

0.36 

11.85 

+0.612 

46.8xl0 6 

FC46 

0.527 

10.4 

+0.684 

46.9xl0 6 

FC70 

0.97 

4.37 

+0.310 

88.77xl0 6 

FC25 

0.242 

19.84 

0.725 

32.22xl0 6 

FC50 

0.434 

13.56 

+5.31 

39.41xl0 6 

FC51 

0.441 

12.89 

-4.58 

38.95xl0 6 


Table 5. Associated Engine Parameters* for these Flight Conditions. 


Flight 

Condition 

Free 

Stream 

Altitude, 

ft 

Free 

Stream 

Mach 

Number 

Inlet Duct 
Exit Static 
Temperature, 
degs R 

Inlet 

Duct 

Exit 

Static 

Pressure, 

psia 

Inlet 

Duct 

Exit 

Velocity, 

ft/sec 

Inlet 

Duct 

Exit 

Mach 

Number 

Mixing Plane 
Total 

Temperature, 
degs R 

Mixing 

Plane 

Total 

Pressure, 

psia 

FC7 

5000 

0.304 

498 

11 

379.6 

0.347 

1050 

23 

FC19 

10000 

0.36 

485.8 

10.2 

345.8 

0.32 

1050 

21.5 

FC46 

24000 

0.527 

443.6 

5.85 

404.3 

0.39 

1045 

14.8 

FC70 

22300 

0.97 

519 

10.65 

464.7 

0.416 

1200 

30 

FC25 

10000 

0.242 

470.1 

8.72 

474.8 

0.447 

1209 

26.3 

FC50 

24000 

0.434 

440 

5.16 

483.3 

0.47 

1154 

16.95 

FC51 

24000 

0.441 

431.8 

5.19 

468.6 

0.46 

1146 

16.74 


&r*The numbers in this table do not represent any particular engine. * ( Formatted. Bullets and Numbering j 
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Concluding Remarks 


A. Aircraft companies, government laboratories, and universities have access to database sets that may include 
flight data, but many of these sets are considered proprietary and are only used by them to assess and improve 
their codes. They hold this information in a closed environment so as to maintain a competitive advantage. The 
CAWAP dataset, being open, allowed researchers in these organizations to predict flight data in a cooperative 
manner under the RTO “umbrella” and facilitate doctoral studies. 

B. There is interest in the international aeronautical community (airframe companies, government laboratories and 
universities) in being able to predict the flow physics measured on a fighter aircraft. 

C. An international team of experts can be assembled with enough patience and institutional support when there is 
a focused common problem of mutual interest that can provide positive payoff for each organization. 

D. The efforts expended by the participating organizations/researchers described here have led to the development 
of improved use or best practices for their flow solvers for the F-16XL aircraft. Moreover, an improvement in 
the ability (Technology Readiness Level) of organizations/researchers to predict complete fighter aircraft flow 
physics at real flight conditions has occurred, in part, due to their participating in this shared and open 
environment. 

E. The preceding leads to the observation that focused datasets should be collected on aircraft configuration(s) of 
international aeronautical community interest without geometrical restriction under the RTO umbrella for 
similar CFD solver improvement. 

Summary 

This paper has traced the F-16XL-1 aircraft and the flight flow-physics data from a NASA-only activity to 
encompass others in the NATO community interested in predicting these data. The Cranked Arrow Wing 
Aerodynamics Project (CAWAP) has been internationalized under the auspices of the scientific arm of NATO and 
the Technology Readiness Level of computational tools has been increased. Sample results obtained by CAWAPI 
facet members have been highlighted to show the breadth of work to be presented by them in their own papers. 
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